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1. Introduction 

In several branches of classical and quantum field theory, as well as in the current attempts to 
develop a quantum theory of the universe and of gravitational interactions, it remains very useful 
to describe physical phenomena in terms of differential equations for the variables of the theory, 
supplemented by boundary conditions for the solutions of such equations [1, 2], For example, 
the problems of electrostatics, the analysis of waveguides, the theory of vibrating membranes, the 
Casimir effect, van der Waals forces, and the problem of how the universe could evolve from an 
initial state, all need a careful assignment of boundary conditions [1-3]. In the latter case, if one 
follows a path-integral approach [4], one faces two formidable tasks: 

(i) the specification of the geometries occurring in the “sum over histories” and matching the 
assigned boundary data; 

(ii) the choice of boundary conditions on metric perturbations which may lead to the evaluation 
of the one-loop semiclassical approximation. 

Indeed, whilst the full path integral for quantum gravity is a fascinating idea but remains a formal 
tool, the one-loop calculation may be put on solid ground, and appears particularly interesting 
because it yields the first quantum corrections to the underlying classical theory (despite the 
well known lack of perturbative renormalizability of quantum gravity based on Einstein’s theory). 
Within this framework, it is of crucial importance to evaluate the one-loop divergences of the theory 
under consideration. Further to this, the task of the theoretical physicist is to understand the deeper 
general structure of such divergences. For this purpose, one has to pay attention to all geometric 
invariants of the problem, in a way made clear by a branch of mathematics known as invariance 
theory [5]. In our paper, however, we shall be concerned with a physics-oriented presentation, and 
hence we here limit ourselves to say that the one-loop divergence receives contributions from the 
Riemann curvature of the background, the extrinsic curvature of the boundary, a gauge curvature 
(if a gauge field is studied), jointly with some matrices occurring in the differential operator and 
in the boundary conditions of the problem. All these ingredients will become sufficiently clear in 
the following sections. 

The recent investigations in Euclidean quantum gravity [2, 6-9] and quantized gauge fields 
[ 10 , 11 ] on manifolds with boundary have shown that boundary conditions involving both normal 
and tangential derivatives of a field arise naturally if one requires invariance of the whole set of 
boundary conditions under infinitesimal diffeomorphisms [2, 6-9] or Becchi-Rouet-Stora-Tyutin 
(BRST) transformations [11]. Thus, shortly after the correct formulae for heat-kernel asymptotics 
with a mixture of Dirichlet and Robin boundary conditions were obtained for the first time [5, 
12 ] after several decades of dedicated efforts by many authors, it became clear that, even just for 
Laplace type operators on manifolds with boundary, the most difficult piece of work in spectral 
geometry [5] is yet to be done. We are thus interested in a framework whose description is as 
follows. 

A smooth manifold M of dimension m is given, with Riemannian (i.e. positive-definite) metric 
g. Moreover, a vector bundle V over M is given with connection V and curvature fl, say. There is 
also a smooth (m — 1)-dimensional boundary dM with the induced metric 7 and the resctriction 
V of V to dM. An operator of Laplace type is a second-order elliptic operator on M which can 
be put in the form (summation over repeated indices is understood) 


P = -g ab V a V b - E 


(l.l) 
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where E may be viewed as a “matrix”-valued function. In general, one studies tensor or spinor 
fields, for which the operator P carries a finite number of discrete indices. The general framework 
for the investigation of one-loop ultraviolet divergences in quantum field theory, consists of the 
heat-kernel approach with the corresponding ^-function regularization [2, 13]. The heat kernel is, 
by definition, a solution of the heat equation 

+ p) [/(*,*';*) = 0 (1.2) 


for t > 0, obeying the initial condition 

lim / dvol(x')U(x, x'\ t)ip(x') = (p(x) (1.3) 

t ^° J M 

with dvol = dx \/detg , jointly with suitable boundary conditions 

[BU(x,x'-,t)] dM = 0. (1.4) 

The functional trace of the heat kernel is obtained by integrating the trace of the heat-kernel 
diagonal U(x,x-,t ) over M, and reads 

Tt L 2 e~ tp = / dvol(x)Tr[/(x,x;f) (1.5) 

JM 


where Tr is the matrix trace. 

In our paper, following [5, 14], we are interested in a slight generalization of eq. (1.5), where 
e~ tp is “weighted” with a smooth scalar function / on M. More precisely, we will be interested in 
the asymptotic expansion as t —» 0 + of the functional trace 

Tr L 2 ^/e _tp j = J dvol(x)f(x)TvU(x,x-,t). 

The results for the original problem are eventually recovered by setting / = 1, but it is crucial 
to keep / arbitrary throughout the whole set of calculations, as will be clear from the following 
sections. 

The asymptotic expansion we are interested in has the form [5, 10, 14, 15] 

OO 

Tr i2 (/e- tp ) ~ (■ i7rt)- m/2 J2t n/2A n/2(f,P,B). (1.6) 

n =0 

The coefficients A n / 2 (f,P,B) consist of two different parts, the interior part C n / 2 {f, P) and the 
boundary part B n / 2 (f, P, £>), i.e. 

A n/2 (/, P, B ) = C n/2 (/, P) + B n/2 (f , P, B). 

The interior parts C n / 2 (f, P) are obtained by integrating some geometric invariants over M and 
do not depend on the boundary conditions. They are calculated up to C 4 in [16]. By contrast, 
the boundary parts B n / 2 (f,P,B) are obtained by integrating some geometric invariants over the 
boundary dM and do crucially depend on the boundary conditions. It is also well known that, for 
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odd values of n, the interior parts C n / 2 (f, P) vanish. In this paper we are mainly interested in the 
boundary parts. 

The first few coefficients for Dirichlet and Robin boundary conditions with 

B d = 1 B r = V n + S (1.7) 

where Vjv is the normal derivative and S' is a matrix-valued function on the boundary, are listed 
in appendix A. The crucial point is that the numerical coefficients multiplying all possible geomet¬ 
ric invariants contributing to the heat-kernel coefficients with the given boundary conditions are 
universal constants , independent of the dimension m. 

As shown in [10], one can consider for the operator P not only the Dirichlet or Robin boundary 
conditions, but also the more general case when the boundary operator takes the form 

b a = + \ (r'v* + wr) + s (1.8) 

where P (i = 1,2,..., m — 1) and S are some matrices which depend on the local coordinates on 
the boundary. Symmetry of the operator P is ensured if P = —P and = S. The consideration 
of boundary operators of the kind (1.8) is suggested by self-adjointness theory [9, 10], string theory 
[17, 18], and BRST invariance [11]. It should be stressed, however, that the results derived in [10] 
only hold when the matrices P commute with each other, i.e. [P,P] = 0 (Abelian case). This 
non-trivial point was not realized by the authors of [ 10 ], and its proof is part of recent work by us 

[19]- 

Our paper, relying on [5, 10, 14], studies the heat-kernel asymptotics for the operator P 
resulting from the boundary operator (1.8). For this purpose, section 2 describes the problems 
occurring in the analysis of Euclidean quantum gravity. The first original results are presented 
in section 3, devoted to the investigation of the totally flat case. Section 4 describes conformal- 
variation methods in a way accessible to physics-oriented readers, with application to the analysis 
of the Ai coefficient in the Abelian case. The structure of A 3 / 2 (f,P,B) and A 2 (f,P,B ) is presented 
in section 5 under three assumptions, the main one requiring again that the matrices P commute 
with each other. Section 6 is devoted to an original reduction technique for heat-kernel coefficients. 
Results and open problems are described in section 7, and relevant details are given in the appendix. 


2. Boundary operator in Euclidean quantum gravity 

As shown in detail in [2, 6, 7, 9, 11], generalized boundary conditions involving also tangential 
derivatives occur naturally in Euclidean quantum gravity. The metric perturbations are symmetric 
rank-two tensors. More precisely, in geometric language, they are smooth sections of the vector 
bundle of such tensors. This bundle has a connection u> e: ab cd = ~^B^ c e ^ a d d \) and a curvature 
Bt e f,ab cd = —2 R e f( a < ' c ^b) d \ and its metric is defined by the equation 

H ab cd = g Ac g d)b _ }_ g abgCd' ( 2 . 1 ) 

The corresponding operator in the covariant de Donder type minimal gauge is then of Laplace type 
(1.1), with the “potential” term of the form 

E ab cd = -2 R {a c b) d + 2R { J c 6 b) V - g ab R cd - R ab g cd + g ab g cd R- 

v ' v ’ m — 2 m — 2 
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Hereafter, N a is the unit normal to the boundary, e® is a local basis for the tangent bundle of 
the boundary, e, is the dual basis on the cotangent bundle, g ab e l a e : > b = 7®-', 7* J being the induced 
metric on the boundary, q ab = g ab — N a N b = e“e^ 7 * J ; Vi = e“V a , and Vat = JV“V„ denotes the 
normal derivative. One has = 0, VjiV = Va/A = Kjej, where K, t j is the second fundamental 

form of a surface dM(r ) obtained by a normal geodesic shifting of the boundary to a distance r 
[5]. Further, we have Vje. t = KijN + tu^-efc, where is defined by this equation. The connection 

V of dM(r ) is defined by = Vj-e, — u^efc and is compatible with the metric 7 , Vi'fjk = 0, if 
the quantities are the Christoffel symbols computed for the metric 7 ij . Here, as usual, we have 

included in the tangential derivatives V both the standard Levi-Civita connection of the boundary 
and the restriction of the bundle connection to the boundary. In the coordinate basis on the 
boundary one has just N a = N a = (1,0,..., 0) and e l a = (0, <5*). A form of the boundary operator 
which leads to boundary conditions completely invariant under infinitesimal diffeomorphisms of 
metric perturbations reads [9] 

b = fiU + (i-u) (hv n + F i Vi + ( 2 . 2 ) 

where /j is an auxiliary dimensional parameter introduced to preserve dimensions, H denotes the 
matrix H = ( H a b cd ), n is a projector of the form 

Uab Cd = q C (a q d b ) (2.3) 

and r® and S are matrices acting on symmetric second-rank tensors, of the form [9] 

r ab cd = —N a N b e i ( c N d ' > + N {a e^) N c N d (2.4) 

S ab cd = —N a N b N c N d K + 2N {a e i b) e^N d ^ [K^ + ^K ]. (2.5) 

As is shown in [9], the operator P with the boundary operator (2.2) is symmetric with respect to 
the metric ( 2 . 1 ). 

We see that the matrices T® (which are, more precisely, endomorphism-valued vector fields on 
the boundary) are antisymmetric and the matrix S is symmetric. Moreover, one finds 

rr j = --Y j N a N b N c N d - N^Pb^^N^. (2.6) 

Therefore, we find that the matrices T® do not commute with each other, 

[r, r j ] = + N {a e> b) e* c N d \ (2.7) 

Further, we compute 

r 2 = -\{m- l)N a N b N c N d - N { J c b) N d \ (2.8) 

This implies that the matrix T 2 commutes with the matrix S 

T 2 S = ST 2 = -(to - l)N a N b N c N d K 
- b) e^ c N^ [K i: j + JijK\ 


(2.9) 
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but does not commute with the matrices T*. Indeed, one has 

r l T j T k = ^N a N b Y j e k{c N d) - (2.10) 

and, hence, 

[r 2 ,r l ] = \{m- 2) [N a N b e^ c N d ^ + IV (o e < 6) IV c iV d ]. (2.11) 

At a deeper level, as shown in [19, 20], these properties reflect the lack of strong ellipticity in 
Euclidean quantum gravity, when the boundary operator (2.2) is used. Without giving any technical 
detail, for which we can only refer the reader to the work in [19, 20], we can say that this means 
that the trace of the heat-kernel diagonal is non-integrable near the boundary. This property 
implies, in turn, that the one-loop calculation is ill-defined, which is much worst than being just 
very difficult. For the time being, one has thus to resort to the analysis of very simple cases, which 
is what we begin in section 3. 

3. Heat-kernel asymptotics in a particular case 

When a general structure is investigated, it is appropriate to show in detail what it implies in 
the simplest possible case. This first step may be eventually used to solve a more general problem. 
This can help to determine at least some class of the terms in the heat-kernel asymptotics. In this 
section we evaluate all terms in all the heat-kernel coefficients that depend only on E, S and T* 
(but not on their derivatives). It is clear that, to calculate these terms, it is enough to replace 
the manifold M by a flat manifold Mo = R+ x dMo, with a flat boundary 8Mq = R m _i, and 
to assume that all the matrices E, S and T' are constant. Note that the manifold M$ is non¬ 
compact and, therefore, the trace of the heat kernel does not exist. However, for a function / with 
compact support, the trace Tr L 2 (fe~ tp ) does exist and reproduces correctly the local results. We 
are actually interested in the boundary contributions, which are effectively present only in a narrow 
strip near the boundary. The difference with respect to the exact result is always asymptotically 
small as t —> 0. In other words, we assume that R a bcd = = K,j = 0, VE = VS = VP = 0. 

Moreover, we also assume that the matrices W and S commute with E: [E. A] = [E, F 7 ] = 0. The 
heat kernel for this problem can be found explicitly. Denoting the coordinates on the boundary 
by x 1 and the normal distance to the boundary by r, and using the Fourier transform in x l , one 
constructs first the parametrix G( A) and then, by the inverse Laplace transform, the heat kernel 
[19]. Omitting this standard calculation one can present the result in the form 

U(t\x, x') = (47rt)- m / 2 e tE |exp —-^(x —x') 2 + U B (t\x, x')j (3.1) 


Q B (t\x,x') = 
x 



d(j du> 

7 j-m/2 


exp 


-c 2 - w 2 + iQ 


+ i(u> + is) 


(r + r') 

~VT. 




where 
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where e is a positive infinitesimal constant. Herefrom, by taking the limit x = x' we obtain an 
integral representation of the heat-kernel diagonal 

U(t\x, x) = ( 47 rt ) -m / 2 e tE [1 + w(t, r )] (3.2) 


where 


w(t, r) = Qs{t\x, x) 



d(j dcu 

yj-m/2 


exp 



UJ 


-|- 2i(co T is) 



1 - 2(r j Q - iy/ts) \lo + is + T j Cj 



(3.3) 


Note that this result is valid in the non-Abelian case when the matrices T 7 do not commute with 
each other and with the matrix S. The only condition is that both T - 7 and S should commute with 
E. 

Herefrom, by integrating the heat-kernel diagonal with a weight function / over Mq we obtain 
the functional trace of the heat kernel 

Ti '£2 [fe~ tp ] = (4nt)~ m / 2 [ e tE Tv (/) 

J M 0 

n poo 

+ (Ant)~ m / 2 Tr e tE / dr w(t,r)f(r). (3.4) 

JdM 0 Jo 

By taking the limit t —> 0 we recover all terms in heat-kernel asymptotics which contain only 
the matrices E, S and T® (without derivatives). It is clear that the first part, in form of a volume 
integral, gives only the interior contributions and does not depend on S and IN , whereas the second 
one gives a purely boundary contribution. Thus, by rescaling r —*• ry/t, jointly with the Taylor 
expansion of f(ry/t) and the integrand in t we find the boundary part of heat-kernel asymptotics. 

The first non-trivial coefficient reads [19, 20] 


A 1/2 (f,P,B)= f Tr/^{1 

JdM 1 

+2 / ,--'„,„ e ^ |e rr '^ - 1 


'R" 


-l 7r( m 


This result is valid in the general case of non-commuting matrices T- 7 . However, the non-com- 
mutativity of the matrices T'' makes the analysis of heat-kernel asymptotics extremely difficult. 
Therefore, we consider hereafter the purely Abelian case when the matrices T J commute with each 
other and also with the matrix S: [T®, T- 7 ] = [T 7 , S] =0. Then the Gaussian integral in the formula 
for A \/2 can be easily evaluated (see also eq. (4.15) in the following section), 


,P,B)= [ Tr/ ^{ 2 (l + r 2 )- 1 / 2 -l}. 

JdM 1 1 J 

This coincides with the result of [10]. We would like to stress that this is valid only when all 
the matrices T J commute with each other, but is not valid under the only assumption that T 2 
commutes with T* (for details, see [19, 20]). 
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Further, in the completely Abelian case (only this case will be considered from now on) one 
w(t,r) =e~ r2/t jl + ^ dz(2 < Sv / t-^r 2 )e[^ 5 “^ (1+r2) “^] | (3.5) 


(which coincides with [10]) and, therefore, 

r oo 00 

/ dr w(t , r)f(r) ~ ]T t (fc+1)/2 ]T Pn,kS n f^ k ~ n) 


(3.6) 


k =0 


0<n<fc 


where = [£,f] r ._ 0 is the n-th normal derivative of /, and p n are some universal functions 
defined by 

1 


Po,k{ u ) = 


k\ 


i r f k + l 


and, for n > 1, 


Pn,k{v) = 


1 


n\(k — n)\ L 


- r-ai )fe (z/) 


2na n -i t k- n (v) r cTn+i ,k—n(y) 


with v = (1 + r 2 )/4. Here a n ^ are the universal integrals (with k,n> 0) 


&n,k 


poo poo 

(v)= / dy dze- y2 - yz ~ uz2 z n y k . 

Jo Jo 


(3.7) 

(3.8) 

(3.9) 


'0 jo 

After making the substitution y = zu, these integrals take the form 


«n,fc(^) = -r 


1 ^ / n + k 


+ 1 


du 


u 


k 


r o 1 . 1 (n+fc)/2+1 

0 [w 2 + u + vy- ” 


(3.10) 


Interestingly, all the functions a n ^ can be expressed in terms of elementary functions. The first 2 
universal functions can be computed explicitly, 


ao.oW = 


1 


2\/l — 4z/ 


Art anh (Vl — 4zz) 


ai,o(^) = 


7T 


2 \/iz(l + c 2\J~u) 


where Artanhz = log The functions a n>k satisfy obvious recursion relations 


Q^OL n ,kiy) — (%n+2,k iy') • 


(3.11) 

(3.12) 

(3.13) 


Moreover, changing the variables y —> \Jvz and z —> y/\fv we hnd the following symmetry relation: 

a n , k (y) = v {k ~ n)/2 a k , n (v). (3.14) 


Integrating by parts, one can get many other relations among the functions a U)k . In particular, 
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One can thus generate all the universal functions a U) k from cto,o and ay,o only, and one finds 


a 2m , 2l - (-1 r +l g^ K ^ 0 ’ 0 


(3.16) 


«2m,22 + l ~ ( — 1 Y 




— /_1 \m+l l—m—1/ 


<*2m+l,2Z ~ ( — l) m v 


—m—1/2 u 


m+l/2 _ 


(3.17) 


(3.18) 


£*2m+l,2Z+l — ( — l) r 


du m I du l 


^) a °'° 


(3.19) 


We list below some of the universal functions p n ^ explicitly, because they will be used in the 
following section. They are (re-expressing them in terms of T 2 for convenience) 


Po,o = 


0 + T 2 


(3.20) 


P i,i = 


(i + r 2 ) 


(3.21) 


Pq i = Artanh 


P2,2 ~ 


Pi,2 - -pr 1 - 


(i + r 2 )3/2 


7 + r 2 


PH,2 ~ ~ 


T + r 2 -1 i 

f 2 4 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


P3,3 ~ 


3(1 + r 2 ) 2 


i ( i 


P2,3 ~ 


r 2 

i ( i 


Artanh \j— T 2 — -- 


(1 + T 2 ) 


7*0,3 = - 1 + 


7*1,3 = -^2 ;^== Artanh yZ-T 2 - 1 


1 \ 1 


r 2 y 


Artanh v — T 2 


12 V + T 2 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


Note that all the universal functions p n ^ are analytic at T = 0, but have a singularity for T 2 = —1. 
As shown in [19-21], this occurs because the boundary-value problem is strongly elliptic only when 
(1 + T 2 ) < 0. 
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In the course of deriving (3.20)-(3.29), we have used also the following formula for the bound¬ 
ary part of heat-kernel coefficients (cf (1.6)), which results from (3.4), (3.5) and the Taylor expan¬ 
sion of e tE : 


r ( n — 1) i 

~ l — 2 —1 n— 2k — 1 i 

= / Tr E 7T^n-2 fc -i^^/ ( 

JdM ,_ n ■ n A,! 


n—2k—j — l) 


(3.30) 


4. Conformal variations and the coefficient 

We now describe briefly a method for the calculation of heat-kernel asymptotics which turns 
out to be very powerful for the boundary conditions without tangential derivatives [14]. This 
method is based on the behaviour of the heat kernel under conformal rescalings, together with some 
simple exactly solvable cases. In this method, one makes a conformal deformation of the whole 
boundary value problem, i.e. the operator P and the boundary operator B, with a deformation 
parameter e, and one studies the behaviour of the heat-kernel coefficients A n / 2 (f, P,B) under this 
deformation. This is used to determine a set of recurrence relations which reduce the evaluation 
of heat-kernel asymptotics to the solution of a system of algebraic equations for a finite set of 
numerical coefficients. We here present shortly the key elements of this analysis. They are as 
follows. 

Let / be a smooth function on M. Then the conformal deformations of the metric and of the 
inward-pointing normal are [14]: g a b(£) = e 2e ^ dab , N a (e) = e~ £ f N a . The conformal deformation 
of the potential terms E(e) and 5(e), and of T*(e), is chosen in such a way that the operator P 
and the boundary operator transform uniformly, i.e. P(s) = e~ 2£ fp, B(s) = e~ £ ^B. One then 
finds, in particular, that T*(e) = e~ £ f T*. The Christoffel symbols of the Levi-Civita connection of 
M rescale as 

r c ab(e) = r c a6 + e {S a c Vbf + <56 c Vaf - 9 a b g cd V d f) . (4.1) 

Further, for the extrinsic-curvature tensor one finds 

Kab(e) = (V„ N b )(e) = e £f (. K ab - eN a V b f - £g ab V N f) (4.2) 

which implies 

K(e) = g ab (e) K ab (e) = e~ £ f [K - e(m - 1 )V N f]. (4.3) 

As shown in [5, 14], three basic conformal-variation formulae hold: 

^A n/2 (l,e- 2£f P,e- £ fB} = (m-n)A n/2 (f,P,B) (4.4) 

^A n/2 (l,P- £ H,e- £ fB) = A n/2 _ 1 (H, P,B) (4.5) 

^A n/2 (e- 2£f H,e- 2£f P,e- ef B) =0 (4.6) 

where H is some smooth function different from /. The formulae (4.5) and (4.6) are used to 
obtain recurrence relations for the interior coefficients, whilst the equation (4.4) can be used to 
determine recurrence relations which, jointly with equations resulting from eq. (4.6) and from the 
boundary conditions, may determine the boundary contributions [5, 14]. More precisely, these 
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properties are not enough to determine the whole heat-kernel asymptotics. For this purpose, one 
needs some additional information which can be obtained by studying some specific cases, e.g. one 
dimension, manifolds with very high symmetry, product manifolds, and so on. In particular, one 
uses extensively a reduction Lemma along the lines of section 6. Hereafter we assume, as already 
specified in sections 1 and 3, that T* commutes with TL Although very restrictive, this assumption 
has a physical counterpart consisting of complex scalar fields coupled to an electromagnetic field. 
Needless to say, this is only a first step towards more interesting models (e.g. Yang-Mills in curved 
space with boundary [19]). 

Before we complete our calculations, it is appropriate to stress why we have considered the 
conformal-variation method. The two main reasons are as follows. 

(i) The emphasis on conformal-variation techniques has led us to elucidate the structure of the 
A 3 / 2 and A ‘2 coefficients in the Abelian case (see section 5), which is naturally obtained within 
the framework of invariance theory, as shown in [5]. Indeed, the conformal-variation technique 
cannot lead, by itself, to the complete calculation of A 3 / 2 and A 2 . Nevertheless, the results of 
section 3 and of the present section will be shown to lead (in section 5) to encouraging progress 
towards the evaluation of A 3 / 2 - Thus, some new insight is definitely obtained. All this appears 
quite interesting, since the progress of knowledge is normally achieved in small steps. 

(ii) There exists already some evidence, in the recent literature, that this method can lead to some 
advancement of knowledge, because our equations (4.12) and (4.13) (see below) have been used in 
[21] to complete an important investigation of the A\ coefficient. 

Note now that the form (1.8) of the boundary operator suggests immediately that the occur¬ 
rence of T* gives rise to new invariants built from the Riemann curvature tensor, the extrinsic- 
curvature tensor of the boundary, E,S and T* (cf [5]). This simple but crucial observation makes 
it possible to start from equation (Al), jointly with 


B 1/2 (f,P,B) = A 1/2 (f,P,B) = (4tt ) 1/2 [ Tr(7/) 

JdM 


(4.7) 


Bi(f,P,B) 



f{b 0 K + b 2 S)+b 1 f, N + finKij FF) 


(4.8) 


as our ansatz for the first 3 coefficients in heat-kernel asymptotics. Our task is to understand how 
to use this ansatz, conformal-variation methods and yet other properties (see below) to evaluate 
7, bo,b\, b 2 and 07. The resulting algorithm could then be used to derive higher-order heat-kernel 
coefficients [5, 14]. 

First, we use the conformal-variation formula (4.3), jointly with the property 




-fKijVT’ - /,\ l' 2 . 


(4.9) 


The identity (4.4) still holds, since all our unknown coefficients are not affected by conformal 
rescalings (see [14] and below). Thus, when n = 2, eq. (4.4) leads to 


j-A^e-^P^B) 


£=0 


(m-2)A 1 (f,P,B) 


(4.10) 
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where the right-hand side is obtained from eq. (4.8), whilst the left-hand side receives contributions 
from eq. (4.9) and from the conformal-variation formulae in [5, 14]. One thus finds an integral 
over dM where the coefficient of f -n should vanish. This leads to (cf [5, 14]) 

—bo(m — 1) — bi(m — 2) + \^b 2 (m — 2) — (m — 4) — a\T 2 = 0. (4.11) 

Since this algebraic equation should hold for all values of m, one finds the equations 

-bo -h + ^b 2 -l = 0 (4.12) 

bo + 2bi — b 2 + 4 — a\T 2 = 0. (4-13) 

It is now clear that conformal-variation methods are, by themselves, unable to evaluate heat-kernel 
coefficients, since eqs. (4.12) and (4.13) only determine 6 0 and or say, in terms of b] and 6 2 . Thus, 
one needs some additional information. This can be provided by considering some simple exactly 
solvable cases. In particular, using the analysis of section 3 we find that the functions 7 , b\ and b 2 
are determined by the functions p n defined in (3.7)-(3.10): 

7 = 7 tt=Po,o h=Qp 0 p 62 = 6 / 91 , 1 . (4.14) 

Zy/TT 

Thus, under our assumptions, we find from the equations (4.12)-(4.14) and (3.20)-(3.22) the 
complete result (in agreement with [ 10 ], bearing in mind that T 2 < 0 ) 


1 \ 2 ,1 

4 Lvi + r 2 J 

(4.15) 

61 = Artanh (a/— T 2 ) — 3 

(4.16) 

62 = 12 

2 (1 + r 2 ) 

(4.17) 

b o- 6 [ ——^.Artanh ( a/ T 2 ) ( ' 1 + r2 J+ 2 

(4.18) 

<Tl P (1 + r2) ). 

(4.19) 


It should be stressed that 7 , 60 , 61 , 6 2 ,07 are universal functions , not affected by the conformal 
rescalings and independent of the dimension m (although this is clear from the procedure leading 
to eqs. (4.11)-(4.13)). One can easily check that, if T —> 0, 7 , 60 , 61,62 reduce to the constant values 
occurring for Robin boundary conditions. All coefficients are analytic in T 2 , in a neighbourhood of 
T = 0. Equations (4.15)-(4.19) express a deep property, because they tell us that, by virtue of the 
occurrence of tangential derivatives in the boundary conditions, boundary contributions to heat- 
kernel asymptotics in the Abelian case are obtained by integrating over dM linear combinations of 
the geometric invariants of the problem, where all coefficients of the linear combinations are now 
universal functions. 


5. Structure of A 3/2 and A 2 in the Abelian case 
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As in section 4, the properties that we are going to derive hold under the following assumptions: 

(i) The problem is purely Abelian (this is the most important assumption), i.e. the matrices T® 
commute with each other: [T*,r j ] = 0. 

(ii) The matrix T 2 = 7 l; r®T J , which automatically commutes with T J by virtue of (i), commutes 
also with the matrix S: [r 2 ,5] = 0. 

(iii) The matrices T® are covariantly constant with respect to the (induced) connection on the 
boundary: ViT J = 0. 

The latter assumption is made to reduce the number of independent invariants. It is thus easy 
to see that, as a consequence of the antisymmetry and the commutativity of the matrices T®, the 
following property holds: 

(r * 1 • • • r ifc )^ = (—i) fc r* 1 ... pA 

Thus, the products of an even number of matrices T® are symmetric, whereas the products of an odd 
number of T® are antisymmetric matrices; for any symmetric matrix X and for any antisymmetric 
matrix Y one has 

Tra(T)T il • • T^+'X = Tra(T)r l1 • • • Y l2k Y = 0 
where a(T) is an arbitrary function of T 2 . 

As we know, the endomorphisms E and S are symmetric matrix-valued functions. Thus, all 
their covariant derivatives are also symmetric. Hence 

Tr a(T)T il • • • T^ 1 E. k ... t = Tr a(T)T il • • • T iah + 1 S\ k ... l = 0. 

By contrast, the bundle curvature is antisymmetric jointly with all its covariant derivatives. There¬ 
fore 

Tr a(r)r ® 1 ■■■r i 2k n. k ... l = o. 

To obtain the form of A 3 / 2 , one has to consider, further to eq. (A4), all possible contractions of 
the matrices T® with the geometric objects of the form fK 2 , fKS , fVK , /VS, fR, fYl and f-jyK. 
Similarly, the general form of A 2 is obtained by taking all contractions of T® terms with geometric 
contributions of the form: i) fK 3 , fK 2 S , fKS 2 , fRK, fYlK, fEI \, fRS , fYlS, fKVK, fSVK , 
fKVS , /SVS, /Wit, /VVS, fVR, /Vfi, fVE; ii) f. N I< 2 , f. N KS , f iN VK, f. N VS , f, N R, 
f-jyYl, and iii) f-NNK (cf eq. (A5)). By virtue of the properties described above we are led to write 
the following formulae (hereafter, K = K\): 

A 3/2 (f,P,B) = A 3/2 (f,P,B ) + (47T) 1 / 2 -'- [ Tr j/L (Aj :j TT J ) 2 

J dM { 

+ ( 7 3 (h’ijVT') K + <j 4 K ll K l j r l T :i + X 1 K ij T i T j S 

+ //.//v, VI 'I-.' + fi 2 R l Uj rri 


+ f f:\EijYY 1 


(5.1) 
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A 2 (f,P,B) = A 2 (f,P,B) + ^~l Trif\a 5 (K ij T i r j ) 3 + a 6 (%r’P) 2 K 

ooU J dM ^ 

+ a-jKijT^ K ri K rl + a 8 KijT i T j K 2 + agKuK 1 ^ K 
+ a w K u K l 3 K rp rVrT'P + auKuK lr K rj rT3 
+ a 12 (K tj r'T J ) 2 S + a 13 (A,,I 'I'-' .S') K + a^KaK 1 ^S 
+ X-.KijVT'S 2 + A 3 A,,r'.STAS- + X A EK l 3 VT^ + X 5 R lNjN T l T^S 

+ XeR^rr^s + x 7 n iN rs + ^R lNjN K\^v l 

+ HiRiNjN^Ti K r iT r T l + ii 5 R lNjN r^K + ^ 6 R l NlN K i3 PP 
+ n 7 R q iqj K\VT l + ii 8 R 1 njT^K + ^R q iqo K lp Y l T^T l T p 
+ nwR^giKijT^ + + /'/, A, / r'l'' 

+ ji 2 K m rr' + £ 3 j? ii | rJ rT''rT J + a 3 [ 5 ,*, s]r + A 2 .s„r'r-' 

+ piRfuj.jfVr* + p 2 R, Nj N;Nr l r J +b 2 n iN , N r 

+ f, N ^(if,,r’P) 2 + 4 3 A ; ,I T'A + PtKuK’j PF 
+ A 1 A (/ r'i' , .s' + /? 6 p liVjAr rT? + 

+ + f-NN^Kij rp |. (5.2) 

Note that we take here for convenience, as our basic geometric invariants, the components of the 
curvature tensors on M but not on dM: Rjjki , R l iij, RiNjN , ^m and . Further, with our 
notation, A n / 2 (f , P, £>) is given by a formula formally identical to the one in appendix A, where all 
the universal constants for boundary terms are replaced by universal functions which only depend 
on r, and are denoted by the same symbol. Moreover, |<Tj, Aj, m, b t , (3 t . , A;, pj j are new families 

of universal functions. We find it worth stressing that many further invariants would occur in 
eqs. (5.1) and (5.2), but their contributions vanish after taking the matrix traces, for the reasons 
described at the beginning of this section. 

The analysis of the conformal-variation formulae resulting from the consideration of eqs. (5.1) 
and (5.2) provides non-trivial information. In particular, we focus on the A 3 / 2 coefficient, and we 
consider the algebraic equations among universal functions which are obtained when all boundary 
coefficients in eqs. (A2)-(A5) are promoted to the role of universal functions, in agreement with 
our notation for A n / 2 (f,P,B). For this purpose, we now insert eq. (5.1) into eqs. (4.4) and (4.6) 
when n = 3. Bearing in mind, from the analysis of eq. (4.11), that an algebraic equation involving 
m, i.e. A + Pm. = 0 Vm, implies that A and B should vanish separately, we find from here the 


following equations for the universal functions contributing to A 3 / 2 : 

\c Q - 2ci + c 2 - c 6 = 0 (5.3) 

—Cq + 2ci — c 2 + 3cg + pi r 2 = 0 (5A) 

1 1 , 

--c 0 + 2ci - 2 c 3 - c 5 + -c 7 = 0 (5.5) 

c 0 - 2ci - c 2 + 2 c 3 - 2c 4 + 3c 5 — c 7 — (t 3 T 2 - p 2 r 2 = 0 (5.6) 
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-c 7 + C 8 - c 9 = 0 (5.7) 

c 7 - 2c 8 + 3c 9 - XiT 2 = 0 (5.8) 

-4c 5 - 5c 6 + |c 9 - /?i T 2 = 0. (5.9) 


The conformal-variation formulae leading to eqs. (5.3)-(5.9) cannot be presented for length reasons. 
One can say, however, that one has to apply repeatedly the formulae in our section 4 and in the 
appendix of [14], Following the method in [14], some of these equations are obtained by setting to 
zero in eq. (4.4) the coefficients which multiply a basis for the integral invariants in eq. (5.1). 

Moreover, relying on section 3, one can express 10 universal functions in terms of the p U) k 
given by (3.20)-(3.29) (see (5.1) and (A4)): 


192 

Co — -^Po,o 

(5.10) 

192 

c 8 — —J=P2,2 

(5.11) 

192 

c 9 — —f=Pl,2 

(5.12) 

192 

C6 — —f=P0,2 

A 

(5.13) 

di4 = 360ppi d 2 o = 360p 3)3 

(5.14) 

360/9Q, l e 9 = 360p 2i3 ei 0 = 360/?i ;3 e 7 = 360/9 0 ,3- 

(5.15) 


In particular, eqs. (A4), (5.1) and (5.3)-(5.13) imply that we have obtained, in a few steps, 11 
equations for the 18 universal functions contributing to A 3 / 2 . Thus, one is really halfway through 
towards the solution of a non-trivial problem in heat-kernel asymptotics (although the restrictive 
assumptions about commutation of the T* matrices imply that harder problems are in sight). The 
task now remains to use lemmas involving product manifolds [14], or different boundary conditions 
(cf our (6.5)), to complete the evaluation of A 3 / 2 . In the case of the A 2 coefficient, the conformal- 
variation method leads to a smaller percentage of useful algebraic equations, and hence we omit 
here any detailed calculation, since it would not be especially enlightening. 


6. Reduction technique for heat-kernel coefficients 

In the calculation of heat-kernel coefficients, any additional information is helpful. In the 
present section we describe a new method which relates heat-kernel coefficients with different 
boundary conditions, and hence can help to reduce the calculations of some difficult unknown case 
to a well studied problem. 

Suppose that A is a first-order differential operator, with (formal) adjoint Af from which the 
second-order elliptic self-adjoint operators P\ and P 2 are built as 


Pi = A^A and P 2 = AA^ 


(6.1) 
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subject to boundary conditions with boundary operators B i and B 2 , respectively. In the applica¬ 
tions, B 1 will be of zeroth-order, and £>2 of first-order. Let 

B 2 = -Brtf. 


Then the operators P\ and P 2 have the same spectrum except for zero-modes, and for any smooth 
function / on M one finds 


d_ 

dt 



Tr L2 (/e- tP2 ) 


^Tr L 2 (<&(/, A)e~ tPl ^j 


( 6 . 2 ) 


where 

= -{A\[AJ}\ - [A,f]A* + [A\f]A. (6.3) 

In the course of deriving eqs. (6.2) and (6.3), we have used the antisymmetry of [ , ], the cyclic 
property of the trace jointly with the identity 


P 2 e~ tP2 =Ae~ tPl 


(6.4) 


Suppose that it is possible to choose the operator A in such a way that <!>(/, A) reduces to 
a convenient function (this occurs, for example, when A is formally self-adjoint or skew-adjoint). 
Then further identities among heat-kernel coefficients are derived: 


(n-m) A n / 2 (f,Pi,Bi) — A n / 2 (f,P 2 ,B 2 ) = A n/2 _ 1 (<I>, Pi,Bi). 


(6.5) 


Equation (6.5) provides a generalization of lemma 3.2 in [14]. If A is a first-order one-dimensional 
differential operator of the form A = d x — b, where b is a real function on a closed interval, 
then A t = —d x — b, P\ = —d x + b x + b 2 , and P 2 = —d x — b x + 6 2 , with Dirichlet and Robin 
boundary operators: B 1 = 1 for Pi and B 2 = d x + b for P 2 . Moreover, $ becomes just a function 
$ = f xx + 2 / 3 , 6 . Hence eq. (6.5) leads, in particular, to the lemma 3.2 of [14]. 

An important problem in heat-kernel asymptotics is, now, whether a generalization of eq. 
(6.5) exists when the boundary operator includes tangential derivatives. In the affirmative case, 
one could combine the result with the conformal-variation formulae, to compute the unknown 
universal functions in the Abelian case (cf [14]). 


7. Concluding remarks 

The naturally occurring question is what has one learned from the new formulae derived in 
sections 3-6. Indeed, the first non-trivial point is having realized that the analysis of the totally 
flat case (section 3) contributes to the evaluation of a large number of universal functions in the 
curved case. This result is not merely technical, but lies at the heart of further progress in the field 
(see (3.20)-(3.29) and (5.10)-(5.15)). Second, in the Abelian case, the matrix P in the boundary 
operator (1.8) is responsible for a large number of new geometric invariants which contribute 
to heat-kernel asymptotics. Moreover, in the integrands of the boundary contributions to heat- 
kernel coefficients, such invariants are weighted with universal functions, independent of conformal 
rescalings of the background metric. This leads to a purely algebraic contribution of universal 
functions to the conformal-variation formulae. Third, the form of the coefficients A 3 / 2 (f, P, B) and 
A 2 (/, P, B) has been obtained for the first time, under the 3 assumptions listed at the beginning of 
section 5. They have been found to involve 49 new geometric invariants, with respect to the case 
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when the matrix T* does not occur in the boundary operator. Fourth, new algebraic equations 
among universal functions have been obtained (see eqs. (4.12), (4.13) and (5.3)-(5.9)). Fifth, the 
generalized formula (6.5) has been derived. This leads to a relation among heat-kernel coefficients 
with different boundary conditions. 

One has now to complete the calculation of the universal functions of section 5. In the 
general case, one has to take into account the covariant derivatives of the matrices T* with respect 
to the Levi-Civita (and the bundle) connection of the boundary, i.e. terms like K^F 1 ^, in the 
heat-kernel asymptotics. In the case of the A\ coefficient, from dimensional arguments it is clear 
that such terms do not contribute. The only possible term in A\ is /TiT*|j, but this vanishes 
because of the antisymmetry of T* (cf [10], where a generalized DeWitt ansatz for manifolds with 
boundary was used). One should also say that covariant derivatives of T* will occur in the heat- 
kernel asymptotics only polynomially and order-by-order, not in more general functions as the 
r* themselves do. The situation is analogous to the appearance of derivatives of the potential in 
Wigner-Kirkwood expansions [22]. 

Moreover, when the matrices T* do not commute among themselves, only the coefficient Ai/ 2 
is known so far, following the recent analysis in [19]. Unfortunately, in that case there is not even 
a general ansatz for heat-kernel coefficients. The solution of this general problem might result 
from the combined effect of applying a number of different techniques [5, 23], and could require 
the development of a few more lemmas (cf section 6). If these last problems could be solved, one 
could get a deeper geometric understanding of one-loop divergences in quantum field theory on 
manifolds with boundary. The one-loop divergence in four dimensions [2] can indeed be obtained 
by setting eventually / = 1 in the formula for A 2 (f,P,B), once that all the universal functions 
have been evaluated. A further crucial problem lies in the analysis of heat-kernel asymptotics for 
non-minimal operators [5] subject to mixed boundary conditions involving tangential derivatives. 
This step has its physical counterpart in the investigation of quantization schemes in arbitrary 
gauges [24], Its consideration is suggested by the results in [19, 20], where it is proved that a 
Laplace type operator on metric perturbations, jointly with local and gauge-invariant boundary 
conditions, is incompatible with strong ellipticity in one-loop quantum gravity based on Einstein’s 
theory. 

All this seems to add evidence in favour of one-loop quantum gravity having a deep impact 
on quantum held theory and spectral geometry [2, 5-16, 19-21, 23-25]. 
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Appendix A 

First, we describe our notation. Our {a, 6,...,} range from 0 through m — 1 and index a 
local orthonormal frame for the tangent bundle of the Riemannian manifold M, whilst our indices 
{i,j, ...,} range from 1 through m — 1 and index the orthonormal frame for the tangent bundle 
of the boundary of M. The semicolon ; denotes multiple covariant differentiation with respect to 
the Levi-Civita (and bundle) connection of M, whilst a stroke | denotes covariant differentiation 
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tangentially with respect to the Levi-Civita (and the bundle) connection of the boundary. If N 
is the inward-pointing unit normal to dM, the normal derivative of a smooth function / on M is 
denoted by /. n = Vjv/- The curvature of the bundle connection is denoted by Q a b- 

For Dirichlet or Robin boundary conditions, the first heat-kernel coefficients read [14] 

A 0 (f,P,B) = f Tr (/) (Al) 

J M 

A 1/2 (f,P,B) = ( 4tt)* / Tr( 7 /) (A2) 

J dM 


A 1 (f,P,B) = - 


Tr 


IM 


f (a\E + a 2 R) 


+ 


Tr 


* dM 


f {boK + &2<5) + bif-N 


(A3) 


A 3/2 (f,P,B) = (4ir)?-^-{ [ Tr\f(c 0 E + c 1 R + c 2 R i NiN + c 3 K 2 
4o4 [ J dM l V 

+ c 4 K i:j K ij + c 7 SK + c 8 S 2 ) + f, N (c 5 K + c 9 S^j 
+ Cef-NN 


(A4) 


A 2 (f,P,B ) = ^^J M fTr(a 3 E’ a ;a + a 4 RE + a 5 E 2 + a 6 R’ a , a 
+ a 7 R 2 + a$R a i ) R ab + a 9 R a b c d,R abcd + 


+ 



f(d x E, N + d 2 R, N + d 3 K% + d A K ij {ij 


+ d 3 EK + d^RK + d 7 R l isunE + dgRiNjNK 
+ d 9 R l Uj K lj + d w I < 3 + duK t jK lJ K + di 2 K i j K : ' l K l i 
+ di 3 ^iV;i + di^SE + di§SR + di^SR^ nin + d\ 7 SK 2 


+ disSKijK^ + d 19 S 2 K + d 20 S 3 + 


+ f-N i^\E + e 2 R + e 3 R i NiN T 64 K 2 + e^K^jK + e 3 SI\ + eg S 2 ^j 
+ f^NN^e^K + eioS'j + ejf ,a - a N 


(A5) 


The numerical constants {c^, 7 , bi, Ci, di, e^} for Dirichlet and Robin boundary conditions (1.7) are 
calculated, in particular, in [5, 14] (no confusion should arise from the use of the symbol y *- 7 for 
the metric on the boundary in section 2 ). 
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